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QUALITATIVE DETERMINATION OF STRATOSPHERIC AEROSOL CHARACTERISTICS 


Included as supplement to this report are descriptions of our approaches in 
noise analysis and removal, and in relating data points to an altitude above tlx? 
earth's surface. These descriptions will be updated as more is learned but other- 
wise will become a part of the final report. Therefore, we invite criticism by 
the technical monitor's office. 

We are able to report encouraging results this month on the determination of 
attenuation coefficients as a function of altitude— the principal thrust of the 
investigation. The results are still tentative due to our continuing efforts to 
scale correctly. However, because of our progress in scaling, tlW limb brightness 
signature has been inverted and the results are indicative of results to come soon. 
It is clear that aerosol information is contained within the data. After treatment 
for noise removal and for geometric distortions inherent in the sensor data acquisi- 
tion geometry and after intensity calibration, the information can bo extracted. 
Figure i shows the comparison between the brightness model and the measured data 
after calibration and scaling. The model is based on the results reported by 
Elfcrmon (1968) (referenced in a previous report). .The nrxlol has Ixjon converted 
to that signal which would be recorded by the S192 sensor. Figure 2 shows the 
inversion of this measured signature as the profile of attenuation coefficients 
versus altitude. For cornparison we include Elterman's models (1964 and 1968) 
in Figures 3A and 3B which show some of the same characteristics although the 
altitude scale is not the same. Both curves represent the attenuation coeffi- 
cient's as turbidity, the ratio of aerosol coefficients to Rayleigh coef f icients. 

It should be noted that a change in altitude scale will change the aiq litu«J.\s nf 
the profile, further work is presently being concentrated in this area. 
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noise: analysis 


A. Periodic Noi.se 


Given a rejjorted 17 cyclcs/scc {vriodic noise and a F192 scan rat' -1 of 94 
scans/sec we searched for a noise periodicity in the neighborhood of 5.5 scans/ 
cycle. Data used in the initial analysis was frem Pass 47, tape 934527, 
channel 18, scans 101-120. 


Because of the gauss i an noise wo were unable to use the measured data directly. 
We fit the data with two exponential curves, one covering the lower altitude 
and the other covering the higher altitude. r Ilie two curves were smoothly joined 
by weighting one curve decreasingly as the other curve is weighted increasingly; 

e.g. , 

.9 Pj 1 .9 P 2 l .7 r 3 ' 

2 2 2 
.1 P x .2 P 2 .3 P 3 . . . 

(.9 P 1 1 + .1 P 1 ? )( .8 P ? l + .2 ? 2 2 ) 


The exponential curves were then used to represent the data for their 
respeo t i ve scans . 

A twenty point curve was derived by selecting a data point and plotting its 
value frem each of twenty adjacent scans. Several curves were obtained by 
repeating the process for every tenth data point. To emphasize differences 
and to minimize effects of amplitude, the lowest value on each twenty point 
curve was subtracted from tlie respective curves, because the magnitude of 
these resultant curves was constant throughout all curves (through all 
altitudes) it was concluded that the noise was additive. 


We found that the noise from the lower altitudes was erratic. This da* a was 
not considered in the derivation cf tlie noise periodicity . From about point 70 
(assigning 1 to the first useful limb data point) thm* curves appear periodic, 
especially within points 90-130. In this region the noise is definitely periodic 
and all curves are very similar to each other, having a rate of about 4.5 scans/ 
cycle. This is equivalent to a frequency of about 20 cycles/sec. 

This approach was repeated using measured data that had been smoothed. The 
results were unchanged. 

It was assumed that points 90-130 were representative of the periodic noise and 
that in the other regions the periodicity was less apparent due to the gaussinn 
noise or our approximations in smoothing or fitting the measured data. 

To find a representative value for the entire scan the bias curves for each 
of the several points (90-130) on the scan were averaged. The periodic effect 
of the noise was then removed by subtracting the bias value derived for each 
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scan from every [joint on that r.cnn. Because wo chose not to assunc either 
a positive or a negative noise addition and bccrusc a constant value, different 
for cadi curve, (the srrvi 1 lest am? >1 Ltude) was subtracted, no zero reference exists. 
The [jcriodie noise was therefore removed to within on additive constant which 
was determined through the scaling process. 

The following charts are representative of the periodic curves discussed above. 
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SCALING & CALIBRATION 


A. Overlap Calculations 


SENSOR 



FOV 

FIGURE B.l 


snisoR 



The S192 field of view (FOV) is the product 
of solid angle subtended by the sensor and 
the distance to the scene. In the case of 
the ear Ui Y - 435 Km; for the limb of the 
earth Y = 2400 Km. 


To calculate the overlap frr*n pixel to 
pixel wd wish to unfold the cone generated 
by the SI 92 sensor. As slmwn in figure B.2 
the optical line of sight rotates nlxxit an 
axis extended an angle a fran the axis, where 
a = 5°32'. The circumference of the circle 
generated by rotating tlve field of view 
around the axis is 

c = 2ax 


where x = h tan a. The scan line encompasses 
11 6°15 ' of the circle. The length, s, of the 
scan line is then 


s 


116.25° „ 
360 u " 


116.25 

360 


2 n x 
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The genera Ic'd cone can be unloldod 
into a two-dimensional figure as shown 
in figure B.3. Let 0 be the angle that 
arc c subtends and lot B be the angle 
that s subtends: 



B = s/y 


116.25 _ x 
160 " v 


116.25 

360 


2n 


sin a 


= .19563 radians 


Note that 


.000182 1240 

pixel 


pixels 

scan 


.22568 


radians/scan 


The difference between .22568 radians/per scan and .19563 radians per scan is 
the conposite overlap for the scan expressed ih angles. Since there are 1239 
overlaps/scan, overlap per pixel is 


.22568 - .19563 
~ 1239 ~ 


. 0242528 mil 1 i rad i nns 


Expressed in meters this is 

.0242528 milliradians x 2400 km = 58.2 meters 
This overlap on the earth wov. 1 d be 

(.0242528) (435) = 10.5 meters 

This overlap ccnpores with the overlap of 6.65 meters reported in the 
NASA documentation. 
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SCALING AND CALIBRATION 


B. Change of Altitude Increment 



FIGlfRE B. 4 


11k? altitude increment is dependent on 
tlK? position of scan relative to the earth. 
This can ho represented as the angle of 
intersection, 'l, between the scan arc and 
the earth's surface. All altitude incre- 
ments Ax | can bo determined from the 
and known characteristics of tlx? scan arc. 


The arc between two adjacent points subtends an angle of a = 116.25 /1239 pts. 
The arc between the first point and any point i subtends an angle of 
B i = T2J9- '116.25,°. 


We wish to express 0^ in terms of 0^ and fh . 

F’rcm figure B. 5 it is shewn that 0^ = 0j - p. 

- °i - j tsf l ,U6 - 25 °' 

Therefore 

Ax. = s sin 0. 

l l 

vdK?re s is the center to center spacing between points. 
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SCAT .INC AND CALIBRATION 


C. Calibration 

The multiplicative sclo factor, a, needed to calibrate the recorded intensity 
is obtained by fitting the measured curve against a inxlol. 'll to model was derived 
from the attenuation coefficients reported by Elteimm in 19G8 ("UV, Visible, 
and IR Attenuation for Altitudes to 50 km, 1%8", Air Force Cambridge Research 
Laboratories) . Details of tlie derivation of the model is discussed in Appendix C. 

As discussed in Appendix A, the removal of the periodic noise also left an 
undetermined additive constant, b, on the intensity scale. These two constants 
were obtained by solving an overdetermi nod set of linear equations to fitv.1 a 
least squares fit of the measured data to the model. 

Thus, 

yi = a x L t b i = 1 , . — ,30 

where = model intensity 

Xj = measured intensity 

The thirty points were selected from along tlie scan line by incrementing every 
20 points from points 1-81 and every' 2 joints from points 100-148. The errphasis 
on the upper altitudes was because the aerosol contribution approaches insig- 
nificance, therefore deviations of tlie data from tlx? model would also approach 
insignificance, assuming a constant Rayleigh contribution. The lower al- 
titudes must be included to s^xie degree, however, to maintain the sumo general 
curve structure. 
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scai.inc and caliiwation 


D. Sky bet Pointinq Accuracy 

Forcing the measured data curve to fit a model ourv-' was necessary partly 
because of tire inaccuracies in tlic Irousekofcpincj data recorded on the Skylrct 
tape. From NASA document PHO-TR 524 Rev A, page B-l, location of the 
intersection point of the sensor pointing vector and the earth model 
ellipsoid is accurate to within a 3.0 n.mi. tolerance. Extending this 
tolerance from an altitude of 435 Km on the earth to 2400 Km on the limb 
gives a possible error of 


x - 3 n.nii. 


1.85 Km 
n.mi. . 


2400 Km 
" 435' Km 


- 30.6 fan 

For our purposes it is important to locate the field of view to within a 
few kilometers, hence the need to scale the data by what one might reasonably 
expect to find as a function of altitude. 

Several attenpts were made to scale using the Skybet information but all were 
unsatisfactory. Not only was the location of the field of ' icw unknown, but 
also tire spatial increment between data points was unknown ; as follows: 

Let SPU’ be tire sensor pointing vector in FCT' , a coordinate system very' 

— 

close to FCT, and lot: SPU bo the sensor point itrg vector in FCT. SC is the 
spacecraft vector in FCT (Earth Centered True). 

A * SPU’ = SPU 

where A is a 3 x 3 notrix which accounts for tire err • ' in SPU' . 


If 0 is tire angle between tire spacecraft vector and the sensor pointing 
vector emanating from the spacecraft, 


y 
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_ si’ll • - sc 
cos 0 = 

'i 


SPU| • I sc] 


(a-spu * ) • - sc 

| A • STU’ I I SC I 


h = spacecraft altitude 

R = radius of earth 

x = altitude of pointing vector 


= (R + h) sin 0 - R 


x 


X 2 X 1 


(R + h) sin [cos 
(R + h) sin [cos 


( (A _ * W) ; -_SC, . 
| A • SPU‘ | |SC| 

(A - SPU') • - SC 

|A • SPU'-J | SC | 


• sin Icos 


-1 


(A • SPU' ) * - SC 

( —— 

| A * SPU' 2 | | SC | 


)1 


winch shcrws that the difference in altitude is not a simple linear 
relationship with the difference in sensor pointing vectors. The unknown 
transformation matrix A containing pointing errors cannot be removed. 
Calculations shew that constructed A matrices with 1% error on the diagonal 
and very snail off-diagonal errors will, as expected, produce the 30 km 
error tolerance predicted with a non-constant change in altitude between 
adjacent data points. Skybet data were for this reason not us‘xl for 
purposes of field-of-view location. 
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